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van Hove's theory of scattering of probe particles by a 
macroscopic target is generalized so as to relate the differen- 
tial cross section for atomic ejection via stimulated Raman 
transitions to one-particle momentum-time correlations and 
momentum distributions of fD trapped gases. This method 
is well suited to probing the longitudinal momentum distribu- 
tions of fD gases in situ, and examples are given for bosonic 
and fermionic atoms. 
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If an ultracold atomic gas is probed by high-energy 
particle scattering so that the struck boson recoils with 
an energy large compared with typical interactions in the 
target system then measurement of the differential scat- 
tering cross section leads to determination of the dynamic 
structure factor [Q-|3j. There are some situations not 
amenable to such a treatment. For example, it follows 
from the Fermi-Bose mapping theorem || that the dy- 
namic structure factor of a one-dimensional (ID) system 
of hard core bosons (Tonks gas JIH) in a tightly confined 
atom waveguide is identical with that of the corre- 
sponding ideal Fermi gas, although their momentum dis- 
tributions are quite different. Here we introduce another 
approach based on a generalization of the clas- 

sic van Hove theory [jl2j of scattering of probe particles 
by a macroscopic target. We apply this very general ap- 
proach to stimulated Raman transitions between two dif- 
ferent hyperfine levels of a ID trapped gas, for example, 
a Bose-Einstein condensate (BEC). Raman transitions 
have been demonstrated as a mechanism for outcoupling 
of coherent atom beams from a BEC ]l3|-|l5|] , and Bloch 
et al. JT(| have employed radio-wave output coupling to 
measure the spatial coherence of trapped atomic gases. 
Here we propose using Raman outcoupling to determine 
single-atom momentum-time correlations and momen- 
tum distributions of ID trapped gases. The relevant reac- 
tive scattering process in the case of a sodium BEC exper- 
iment is photon + BEC — > photon + boson + BEC 
where photon denotes a photon absorbed from an inci- 
dent laser beam, photon denotes a photon emitted into 
a second overlapping laser beam via a stimulated Ra- 



man transition, boson denotes an ejected atom in the 
untrapped m — magnetic sublevel, which originally 
belonged to an iV-atom BEC of atoms in the m = — 1 
sublevel, and BEC' denotes an (N — l)-atom excited 
state of the residual BEC. The same scheme applies to 
the experiments on Bragg spectroscopy of sodium BECs 
Jl7],[t8| but there the internal state of the scattered atom 
boson remains in the trapped state m = — 1. It is this 
distinction that allows our scheme to yield the momen- 
tum distribution directly by examination of the ejected 
atoms in contrast to Bragg spectroscopy which yields the 
dynamic structure factor. 

Another recent proposal [[l9| for measuring the one- 
particle density matrix via the reactive scattering process 
impurity + BEC — > impurity + boson + BEC simulta- 
neously assumes a high energy incident impurity (Born 
approximation) and an S-wave pseudopotential impurity- 
boson interaction justified only for low energy collisions, 
and requires measurement of a differential cross section 
involving two outgoing particles. Some recent work pc| ] 
on the use of photoionization measurements for deter- 
mination of correlation functions of the Bose field used 
a different scheme (tightly focused laser beam, measure- 
ment of photoelectron energy but not angular distribu- 
tion) with emphasis on temporal rather then spatial cor- 
relations. 

Derivation: A crucial point p"l| ] is that in order to 
probe one-particle correlations one should use a reactive 
scattering process which effectively removes one particle 
from the target system; this leads to the hole propagator 
and a related one-particle momentum-time correlation 
function. For the case of the above-described Raman pro- 
cess, we use the following simple choice of asymptotic ini- 
tial and final states: \i) — | laser) |B,;) = U\Bi) and |/) = 

"l, m =o\ laser) |B/) = 4,m=o^l B />- The laser state vec- 
tor is a single-mode coherent state with respect to each 
of the two Raman beams, one with wave vector ki and 
polarization index Ai and the other with wave vector k2 
and polarization index A2 . Its most convenient represen- 
tation for our purposes is | laser) = Z7|0) where the uni- 
tary operator U = e F with F = c\b\ x + c-^b 



k2>2 



h.c. 
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The mean photon occupation numbers are n,j = \Cj\ 
and the occupation number dispersions are */nJ, negli- 
gible compared to rij ^> 1. For concreteness in nota- 
tion we consider the case of a sodium BEC target but 
stress that the treatment is not restricted to bosonic 
atoms or equilibrium states. Then |Bj) is the iV-boson 



1 



state vector describing N magnetically trapped bosons 
in their electronic ground state and magnetic sublevel 
m = — 1, including all effects of their mutual interac- 



The derivation closely parallels our previous work [[Ll| . 
The exact transition rate from a specified initial state to 
a specified final state is (2ir /Ti)\(f\T\i)\ 2 8(E f — Ei) where 



tions, ajc m=0 creates an atom with wave vector k in Ef — Ei = E-q — E-q — Tko with Tilo = hc(k\ — k 2 ) 



the untrapped m = magnetic sublevel, and |B/) is an 
(N — l)-boson final state of the residual BEC after the 
m = atom has been ejected. The geometry of this pro- 
cess is shown in Figure [I] where we show a cigar shaped 
BEC of length L along the longitudinal direction x. We 
assume that the BEC is tightly confined in the trans- 
verse plane tt = (y, z) using an atom waveguide. At 
resonance the energy difference Tic{k\ — k^) between the 
photon absorbed by the atom (number 1) and that emit- 
ted (number 2) is equal to the energy splitting A between 
a trapped atom with m = — 1 and that of the untrapped 
atom with m = (A includes both the Zeemen split- 
ting and difference in confinement energy between the 
two states). The net momentum h(ki — k2) transferred 
to the atom will be in the forward direction z as shown 
provided that the indicated angles of the two beams sat- 
isfy sin 62/ sin 6\ — ki/k^ ~ 1. The requirement is then 
6 2 = 0i +e with e w [(fci - k 2 )/k 2 ] tan#i < 1. In the de- 
sign of experiments the polarizations Ai and A2 of the two 
laser beams must be chosen appropriately, and a suitable 
detuning must be used to inhibit spontaneous transitions 
]l4j . Inclusion of these details changes the differential 
cross section only by a multiplicative factor and will not 
be considered here. 




FIG. 1. Geometry of the Raman outcoupling 

The appropriate Hamiltonian is H = Hb + Ho P h + 
Hoa + V where Hb is the second quantized many-body 
Hamiltonian of the trapped m = — 1 BEC atoms includ- 
ing their mutual interactions and V is the interaction of 
the TV-atom system with the quantized electromagnetic 
field. The asymptotic (free) Hamiltonians for the inci- 
dent photons and outgoing m = atoms are -f/opft = 
KcfaNi + k 2 N 2 ) and H 0a = EiJ(^ 2 /2m) + A]iV k , m=0 
where Nj = b^.^b^^ and iV k ,m=o = al m=0 a ktm=0 are 
the occupation number operators for photons and m = 
atoms in their electronic ground state and A is the en- 
ergy difference between m = — 1 and m = atoms. The 
coherent laser state | laser) is not an exact eigenstate 
of Hoph , but its energy dispersion is negligible since the 
photon numbers are large. 



{h 2 k 2 /2m) - A, the energy transferred to the BEC; the 
momentum transferred is Ti(\ = ?i(ki — k2 — k). Here T 
is the exact transition operator T — VCt + = Cl~V where 
Q ± are the Moller wave operators U (0, =F°o) and U is 
the adiabatically switched evolution operator The 
angle-energy doubly differential cross section is then iFy] 



p(0,0) / dt{U^f\t)N Km=0 f{Q)U) B (1) 



d 2 a 
dQdiu 



where (--^b denotes an average over a statistical en- 
semble of BEC initial states |Bj) (reducing to an N- 
boson ground state average for temperature T — ► 0) 
and the prime indicates that this expression is to be 
evaluated for an outgoing m = atom with k vec- 
tor in the direction (9, </>) and with magnitude k — 

\J (2m/h 2 )[hc(ki — k 2 ) — A — hu>] consistent with a spec- 
ified energy transfer hut to the BEC. Here p(9, <j>) is the 
density of final states for an observation distance D from 
the ID gas which is much larger than the length L of 
the trapped gas, (j) is the angle of atomic ejection in the 
transverse plane, and 6 the angle of ejection measured 
out of the transverse plane (see Figure |l|). T^(t) is the 
Heisenberg operator 

fUt) = e lt (HB + H 0ph +H 0a )/hfj e -it{H B + H 0ph +H 0a )/h 

and the BEC final states have been eliminated by an ex- 
act closure summation J^B 1^/) (^/ 1 = ^ after introduc- 
ing the integral representation of the energy delta func- 

The terms in T(0) = T contributing to the Raman 
process have the structure 



2^ fo k 2 A 2 a k ^. m=0 ( k 2A 2 ,k 2 |r|k 1 Ai,k 1 )a k ^ 

k 1 k 2 



-l°kiAi 



(3) 

To obtain the relevant terms in T^(t) we take the her- 
mitian conjugate of (3) and propagate each of its anni- 
hilation and creation operators to time t in accordance 
with (2). Heisenberg propagation of b^ kx and a k ' m=0 
gives only trivial phases since the target BEC contains 
only atoms with m = — 1, whereas the propagation of 
a k ' m= _i is nontrivial and leads to inclusion of effects 
of energy and momentum transfer to the target BEC. 
Noting that the T-matrix elements (^2^2, k 2 |T|kiAi, k x ) 
contain a momentum conservation factor S, , , ' , , , / 

ki+k 1 ,k2+k 2 

and performing Wick's theorem algebra after evaluating 
the coherent laser unitary transformation U^ 1 ■ ■ ■ U in 
Eq. (1), one finds 
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d 2 a 
dfldu 



ni{n 2 + l)p(6,<j>) 
x |(k 2 A 2 ,k|T|k 1 A 1 ,k-k 1 +k 2 )| 1 



k-k!+k 2 



(t)a k _ kl+k2 (0)) B e^ t dt (4) 



where the annihilation and creation operators refer to 
atoms in the BEC target, their m = — 1 subscript hav- 
ing been dropped for simplicity in notation. Physically, 
Eq. ([|) reflects that fact that outgoing m = atoms 
with momentum hk. originate from transfer of momen- 
tum 7i(ki — k 2 ) from the Raman photons to the target 
bosons of momentum h(k — ki + k 2 ) and the differential 
cross section at the corresponding angle is related to the 
number of such target bosons. 

ID trapped gases: The expression (||) is exact. We now 
introduce approximations appropriate to tightly confined 
ID gases where only the ground transverse mode is oc- 
cupied 0. When an atom is transferred by the Raman 
transition from the m = — 1 trapped state to the m = 
untrapped state it is ejected and undergoes radial ex- 
pansion in the transverse plane yt = (y, z) in addition to 
the momentum kick ?i(ki — k 2 ) in the forward direction 
z from the Raman process. In the limit that the trans- 
verse mode u 9 (yt) of the trapped atoms has a width wt 
much less than both the length L of the sample and the 
Raman laser wavelength, the ejected atom wave function 
will predominantly propagate out radially as a cylindri- 
cal wave concentrated along the direction 6 = 0. The 
density of final states p(6, may then be replaced by a 
constant (we assume the atoms can be detected before 
gravity significantly changes the final density of states). 
In addition, writing the field operator for the trapped tar- 
get atoms as ip(r, t) = u g (rx)v(x, t) to reflect their single- 
transverse mode structure, and setting k = k^+xfea, with 
k x = ksm.(0) the longitudinal component of the wavevec- 
tor, we find, for example 



ak-k!+k 2 0) = / d A rip(r,t)e 
= u g (K T )c(k x ,t) 



-i(k-ki+k 2 )-r 



(5) 



Here u 9 (Kr) is the Fourier transform over the trans- 
verse plane of the transverse mode evaluated at Ky = 
(k^ — ki + k 2 ), that may be approximated by its value 
%(0) at zero momentum under tight confinement con- 
ditions 2-k/wt » |kx — ki + k 2 |, and c(k x ,t) = 
J dxv(x, t) exp(— ik x x) is the Fourier transform of the 
longitudinal field operator. Putting this together in Eq. 
(f|) and treating the T-matrix elements as constants (see 
below) yields an expression proportional to the Fourier 
transform of the longitudinal one-particle momentum- 
time correlation function: 



d 2 <7 

dQduj 



{c\k x ,t)c(k x ,0))' B 



'dt 



(6) 



where a is a constant. Thus, measurement of the doubly 
differential cross section for the ejected atoms gives access 
to a one-particle momentum-time correlation function of 
the trapped atoms. 

Further simplification is possible if measuring the en- 
ergy transferred to the BEC is not feasible, and only the 
angular distribution is measured. We consider that the 
ID gas is confined longitudinally in a harmonic trap of 
frequency lot, in which case the maximum energy trans- 
fer to the BEC that can occur is hui max = NTllut- If we 
arrange that the Raman transition is detuned such that 
fou max <C Ti8 with (5 = c|ki — k 2 | — A, then for all relevant 
energy transfers k ss fc — \J (2mS/h), and the implicit 
dependence of (•• -) B in Eq. @ on uj may be neglected. 
Then the integral over to yields a delta function of time, 
and one obtains an expression for the singly differential 
angular cross section in terms of the static longitudinal 
momentum distribution of the target: 

^ « (3(c\k o sm(9),0)c(k o sm(6),0)) B , (7) 

where (3 is a constant. Measurement of the differential 
cross section as a function of angle 9 out of the transverse 
plane will then yield the structure of the longitudinal mo- 
mentum distribution of the target atoms. For a longitu- 
dinal trap of frequency u>t the characteristic longitudinal 
momentum is k osc = 2ir/x osc , with x osc = yJfi/mUT the 
ground state width. If ko >> k osc then we may approx- 
imate fco sin(6>) pa k()8 for all angles of interest, in which 
case the differential angular cross section will directly 
mirror the momentum distribution. 

Evaluation of T-matrix: Although the expressions (6) 
and (7) do not require explicit evaluation of the T-matrix, 
it is helpful to examine its leading term to verify that 
it can be absorbed into constant prefactors. The Ra- 
man process is of second order in the interaction be- 
tween the atoms and the laser field, so we approximate 
T by the second term in its Born series: T w H GqH . 
Here H represents the interaction of the atomic elec- 
trons and nuclei with the quantized electromagnetic field 
and G = (E + ie - Hq)- 1 with E = E. L = E f (re- 
membering the energy delta function in the cross section) 
and e — * 0+. Hq is the Hamiltonian of free atoms and 
quantized electromagnetic field. We use the electric field 
gauge with quantized electric field operator 



E(r) 




ekA&kAe 



ik-r 



h.c 



(8) 



where e k A are the unit polarization vectors and V is 
the quantization volume. It is assumed that the Raman 
lasers are tuned close to resonance for a single S-P tran- 
sition from the atomic ground state to an electronically 
excited state v with electronic excitation energy e„o, i-e., 
Tick\ « hcki « e„o- Then the dominant contribution to 
the T-matrix element in Eq. (4) is 



3 



(k 2 A 2 ,k | T^X^k-^ + k 2 ) ^ -(hc/V)^k^ 

x (e k2 A 2 • d 0r /)(e kl Ai ■ d u0 )/ (e u0 - hcki), (9) 

where d„o is the corresponding dipole matrix element. 
The contributions H 2 k 2 /2m (kinetic energy of ejected 
atom) and A (Zeeman splitting) to the energy denom- 
inator have been dropped since they are much less than 
the electronic transition energy. Assuming that the Ra- 
man laser parameters are held constant, the expression 
(^) is indeed a constant. 

Discussion: Figure ^| shows numerical results for the 
angular cross section versus angle for both a ID system of 
N = 10 hard core bosons |^2| in a harmonic trap (dashed 
line) and the corresponding system of non-interacting 
fermions |^3| (solid line), where we chose lot = 2tt x 1 
Hz, S — 2n x 10 3 Hz, so that k = 7.3k osc . 
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FIG. 2. Angular cross section versus angle sin(#) « 9 for 
N — 10. The dashed line is for the ID gas of impenetrable 
bosons and the solid line is for the corresponding system of 
non-interacting fermions 

These results highlight the difference between the mo- 
mentum distributions for Fermi and Bose gases: The cen- 
tral peak characteristic of a Tonks gas gets sharper as N 
increases ]5],^2| , whereas the fermionic momentum distri- 
bution, which shows the rounded Fermi sea due to the 
trapping potential with Friedel oscillations superposed 
p3fl , broadens with increasing TV. We remark again that 
this difference in angular distribution arises even though 
both the fermionic and bosonic atoms have the same 
dynamic structure factor. This highlights the need for 
a measurement technique that accesses the momentum 
distribution directly. One could measure the momentum 
distribution by simply releasing the atoms from the trap 
and letting them expand, but this leaves open the pos- 
sibility of many-body interactions obscuring the results 
as the gas expands. The momentum distribution could 
also in principle be reconstructed from the single-particle 
reduced density matrix px(x,x') as measured in the ex- 
periment of Bloch et al. |n| , and it will be of interest to 
see how our approaches compare. 
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